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: ABSTRACT 

00 ■ We present the results of ground-based imaging spectroscopy of the [Ne ll] 12.8 /im 

O ; hue emitted from the ultracompact (UC) H ll regions; W51d, G45.12+0.13, G35.20-1.74 

^ ' and Monoceros R2, with 2" spatial resolution. We found that the overall distribution of 



O 

^ ' the [Ne ll] emission is generally in good agreement with the radio (5 or 15 GHz) VLA 

p distribution for each source. The Ne"*" abundance ([Ne"^/H+]) distributions are also 

Q ! derived from the [Ne ll] and the radio maps. As for G45.12-I-0.13 and W51d, the Ne"*" 

^ I abundance decreases steeply from the outer part of the map toward the radio peak. On 

^ ■ the other hand, the Ne"*" abundance distributions of G35.20— 1.74 and Mon R2 appear 

rather uniform. These results can be interpreted by the variation of ionizing structures of 
neon, which is primarily determined by the spectral type of the ionizing stars. We have 
^ I evaluated the effective temperature of the ionizing star by comparing the Ne"*" abundance 

d • averaged over the whole observed region with that calculated by H ll region models based 

on recent non-LTE stellar atmosphere models: 39,100 tHo'^ K (07.5V-08V) for W51d, 
37,200 K (08V-08.5V) for G45.12+0.13, 35,000-37,600 K (08V-09V) for 
G35.20-1.74, and < 34,000 K (< BOV) for Mon R2. These effective temperatures are 
consistent with those inferred from the observed Ne"*" abundance distributions. 
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1. INTRODUCTION 

Ultracompact (UC) H II regions, which are known 
to be bright compact Galactic radio sources, are 
thought to be powered by massive young stars stiU 
embedded in their parent molecular clouds (Wood & 
Churchwell 1989 ; hereafter WC). They are very small 
(r < 0.1 pc) and dense (high electron number den- 
sity : Ue > 10* cm^"^), and have a large emission 
measure (EM > 10^ pc cm^^). The spectral type of 
the ionizing star is uncertain for most of the UC H 
II regions because the ionizing stars are heavily ex- 
tinguished by the dust clouds surrounding the UC 
H II regions, which prevent optical and even near- 
infrared observations even though the ionizing stars 
are very luminous. On the other hand, mid-infrared 
fine-structure lines do not suffer so heavily from ex- 
tinction by the dust, and thus are good probes of the 
UC H II regions. 

Since the warm dust which surround the UC H 
II regions re-radiates nearly all of the absorbed stel- 
lar photons in the far-infrared, the UC H II regions 
are known as strong Infrared Astronomical Satellite 
(IRAS) point sources, with flux densities > 10* Jy at 
100 ^m. WC pointed out that the total FIR lumi- 
nosity derived from the IRAS flux is larger than that 
expected from the spectral type of the ionizing star 
estimated from the radio continuum flux for most of 
the UC H II regions they observed. Two possible in- 
terpretations were presented by WC; 1) strong dust 
absorption of ionizing photons inside the UC H II re- 
gion, and 2) presence of a stellar cluster including 
a large number of non- ionizing low-mass stars, which 
substantially contribute to the total luminosity. How- 
ever, due to the uncertainty in the dust extinction in- 
side the UC H II region, it is diflacult to distinguish 
which effect is dominant. 

In this paper, we propose the [Ne II] line at 12.814 
Hm as a probe of the spectral type of the ionizing 
star. The ionization potential of neon (Nc° Ne+ : 
21.56 eV) is higher than that of hydrogen, and there- 
fore, this line is a good tracer of early-type (especially 
BO or late-0) stars which form the UC H II regions. 
Moreover, the [Ne II] flux is not sensitive to Ug as long 
as Tie < Ucrit (= 3.6 X 10^ cm~^; the critical density 
of the line; Lacy, Beck, & Geballc 1982), and arc thus 
good probes of the ionic abundance, even for dense 
UC H II regions. In addition, since the line is within 
the atmospheric window, the A^— band (8—13 /xm), it 
is accessible from ground-based telescopes. 



The sizes of UC H II regions are very small; the 
typical source size of ~ 0.1 pc corresponds to 2" on 
the sky if the source is located at 10 kpc away from 
us. Thus the [Ne II] line observations in the previ- 
ous work could not spatially resolve the source due to 
their large beams (7—30"). The 2" spatial resolution 
of the present work enables us to obtain the spatially 
resolved [Ne II] images of the UC H II regions, which 
can be directly compared with the radio VLA images 
(WC; Kurtz, Churchwefl & Wood 1994). 

In the recent work by Watarai et al. (1998; here- 
after WMTM), we presented [Ne II] observational 
resuhs for two UC H II regions, G29.96-0.02 and 
W51d, using the same instrumentation, and esti- 
mated the effective temperature of the ionizing star 
from the Ne"*" abundance, by using H II region mod- 
els based on recent stellar atmosphere models. Here, 
we present new observational results for four UC H II 
regions including W51d, for which observations were 
performed over a wider area than that of WMTM. 
Because of the increased number of observed sources 
in this work, we can discuss the relation between the 
spectral type of the ionizing star and the Ne"*" distri- 
bution of the UC H II regions. 

The selected four UC H II regions are W51d, 
G45. 12-^0.13, G35.20-1.74, and Mon R2 from the 
WC source list. W51d, also known as W51 IRS 2, 
is the brightest region in the W51 star forming re- 
gion. The total FIR luminosity is estimated to be 3 
X 10® Lq (Erickson & Tokunaga 1980), which is a 
factor of 3 brighter than the bolometric luminosity of 
an 04 zero-age main sequence (ZAMS) star, whereas 
the Lyman continuum luminosity {N^yc) derived from 
the radio flux corresponds to that of an 05.5 star. 
G45. 12-1-0. 13 is a well-studied UC H II region at many 
wavelengths. Lumsden & Puxley (1996) showed very 
good positional coincidence between the Br7 peak and 
the radio continuum peak. WC noted that an 04 (sin- 
gle) star or 05 (cluster) star is required to explain the 
FIR luminosity (1.7 x 10® Lq), whereas a single 05 
star can account for the Lyman continuum luminosity 
(log NLyc 49.53). G35.20-1.74 is part of the W48 
complex. According to WC, the total FIR luminosity 
(2.8 X 10'"' Lq) measured by IRAS is larger than that 
expected from the spectral type (07.5) derived from 
the radio observation. As for Mon R2, its FIR lumi- 
nosity is dominated by IRS 1. Downes ct al. (1975) 
estimated the spectral type of the ionizing star to be 
BO ZAMS or earlier star based on the radio observa- 
tion. Howard, Pipher, & Forrest (1994) also reported 
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that the ionizing star located at IRS Isw (see Figure 
1 or Table 4 in their paper) is BO— B3. The Galacto- 
centric distances (Rg) of these sources are similar to 
that of the Sun {Ra = 6—9 kpc), and hence we as- 
sume a constant elemental abundance of neon equal 
to the local interstellar abundance (hereafter LIA) : 
9.9 X 10-5 (Simpson et al. 1998). 

This paper is organized as follows. In §2, we de- 
scribe the observations, and the data reduction. In 
§3, we present the [Ne II] images, and derive the Ne"*" 
abundance map. In §4, we discuss the effective tem- 
perature for each source by using calculations based 
on the latest stellar atmosphere models in relation to 
the ionization structure. We also give some remarks 
for individual sources. Finally, §5 gives the summary 
of this work. 

2. OBSERVATIONS 

The [Ne 11] 12.8 fim line images were obtained 
at the Wyoming Infrared Observatory (WIRO) us- 
ing a liquid helium cooled Mid-Infrared Fabry-Pcrot 
Imager (MIRFI; Watarai et al. 1996) developed at 
Nagoya University. MIRFI incorporates a Si:P 5 x 
5 detector array sensitive to mid-infrared wavelength 
(< 28 /Ltm), and a Fabry-Perot interferometer with 
spectral resolving power A/AA = 1600 (190 km in 
velocity). The plate scale with the WIRO 2.3-m tele- 
scope is 1'.'7 pixel" ^ and so the field-of-view (FOV) is 
10'.'2 X 10'.'2 including gaps between pixels. 

The observations of the UC H II regions were car- 
ried out on 4— 11 October 1997. The telescope beam 
was mainly pointed at the radio peaks. However, 
since the FOV of MIRFI is not large enough to cover 
the whole H II region, we carried out a raster map- 
ping. For each source, the raster mapping was carried 
out at various spatial positions around the radio con- 
tinuum peaks. One frame of data was taken at each 
raster position. Each raster position was accurately 
determined by position analysis using observations of 
the guide stars. The spatial sampling interval of the 
raster frames is 0'.'3 — O'.'S. After the flat-fielding using 
the blank sky data, the final images were made up by 
mosaicing of the images taken at all raster positions. 
Finally we could obtain wider spatial coverage images 
of each source. Details of observations arc shown in 
Table 1 and 2. The secondary chopping frequency 
is 2.5 Hz, and the chopper throw was 40" along the 
east-west direction. Furthermore, nearly all data were 
taken by nodding of the telescope beam by 40" along 



the cast- west direction in order to remove the residual 
sky emission features which were visible in the data 
with a single chopping beam. 

The [Ne II] spectral data were acquired by wave- 
length scanning of Fabry-Perot interferometer in 
"LINE SCAN" mode (Watarai et al. 1996). In this 
mode, the interferometer scans around the [Ne II] line 
wavelength in 16 steps over the wavelength range of 
0.034 /im (770 km in velocity) with a step size of 
0.0021 nm (48 km s-^). It takes 12.8 sec to perform 
a LINE SCAN. We define 40 LINE SCANs, half of 
which are taken at each nodding position, as 1 frame, 
hence the total observing time of 1 frame is 512 sec. 
The 3a detection limit was estimated to be 6 x 10"^° 
W cm"^ pixel"^ for 1 frame observing time. 

The absolute flux scale was calibrated using the 
continuum fluxes emitted from the mid-infrared 
bright stars a Ori, x Cyg, and a Tau. The flux den- 
sities at 12 /im of these stars are assumed to be 2820 
Jy, 850 Jy and 450 Jy, respectively. The flux of a Ori 
is taken from the IRAS Low Resolution Spectrometer 
(LRS) data (Olnon et al. 1986), the flux of x Cyg is 
taken from Gillett, Low, & Stein (1968), and the flux 
of a Tau is given in the mid-infrared spectra mea- 
sured by Cohen & Davies (1995). The uncertainty of 
the absolute flux calibration is estimated to be ± 12 
% including the uncertainty in the line-fitting proce- 
dure. 

The observed spectra at the [Ne II] peak positions 
are shown in Figure 1. Each spectrum has been fitted 
with a constant continuum and a Lorentzian line pro- 
file. The spectral resolution of MIRFI is larger than 
the intrinsic [Ne II] line width expected for UC H II 
regions (< 100 km s"^), and so the derived line pro- 
files represent instrumental ones. The central wave- 
length and the wavelength scale were determined by 
the measurement of two nearby bright atmospheric 
emission features (12.832 /im and 12.870 /;,m). As a 
result, the FWHM line width and the central wave- 
length were determined to an accuracy of Av < 20 
km s"^ (3cr). 

MIRFI is equipped with an iJ— band monitor for 
guiding the source into the MIR detector FOV and for 
checking the pointing accuracy. During each observa- 
tion, bright SAO stars near the sources were pointed 
before and after every [Ne II] line observation of the 
sources. The resulting uncertainty in the position de- 
termination was estimated to be < 1". For the posi- 
tional data of the guide stars, we used the latest HST 
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guide star catalog.^ 

3. RESULTS 

3.1. [Ne II] Images 

As we described before, the raster mapping was 
carried out for each source. The hnc fluxes at same 
sky positions between the frames taken at different 
raster positions agree well within the statistical er- 
rors. Therefore, the relative flux calibration among 
all raster positions is reliable. The number of raster 
positions (= number of the frames) obtained for each 
source are listed in Table 2. Mon R2 was observed 
over a lot of frames, but only one-third of the radio 
emitting region was covered, since Mon R2 extends 
to more than 25" wide. Furthermore, we took frames 
at only 5 raster positions for G35.20— 1.74, hence the 
whole emitting region could not be covered. 

The final [Ne II] images smoothed by a 1'.'7 FWHM 
(1 pixel of MIRFI detector array) Gaussian profile and 
corrected for the interstellar extinction are shown in 
Figure 2. The [Ne II] peak positions of all sources ex- 
cept for G45. 12-1-0. 13 show good agreement with the 
5 or 15 GHz radio continuum peaks within the posi- 
tional error. As for the correction for interstellar ex- 
tinction, the following values are used; Ay = 20 mag 
for W51d based on J, H, K, and radio observations by 
Goldader & Wynn-WiUiams (1994), Ay 24 mag for 
G45.12+G.13 from the Br7 imaging and the radio con- 
tinuum observations by Lumsden & Puxley (1996), 
Av - 23 mag for G35. 20-1.74 derived from CO ob- 
servations (Vallee & MacLeod 1990) and from the 
Bra/Br7 ratio (Woodward, Heifer, & Pipher 1985), 
and Ay ^ 23 mag for Mon R2, from the Br7/Bra 
ratio with a relatively large beam (28") by Natta et 
al. (1986). We assume a uniform extinction over the 
whole area of each source. 

However, the [Ne II] line is on the 10 /j,m silicate 
feature and the extinction law from NIR to MIR wave- 
length is still uncertain, especially toward UC H II 
regions. Moreover it may vary with location in the 
Galaxy. Ay also may have spatial distribution inside 
the sources. Hence the uncertainty due to extinction 
should be an additional ~ 20 %. 
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3.2. Physical Properties and Ne+ Abundance 
Maps 

In order to obtain the Ne"*" abundance ([Ne+/H+]), 
we derive the physical properties using the observed 
[Ne II] flux and the 5 or 15 GHz radio continuum 
flux. The procedure is described in WMTM in more 
detail. First, we calculated the brightness tempera- 
ture (Tb) and the optical depth (r) of each source from 
their radio continuum fluxes, and then the emission 
measure (EM) was derived using an expression of the 
free-free radiation including a correction for the opti- 
cal depth effect. We estimated a dependence of the 
electron temperature to EM assuming Case-B recom- 
bination. The EM changes only a few percent under 
the physical condition of typical UC H II regions, even 
if Te changes from 5000 K to 15000 K. Accordingly we 
have considered the error of EM is mainly caused by 
systematic and absolute flux errors. Since the radio 
fluxes are quoted from the same WC survey data, the 
relative errors in the radio fluxes among the sources 
should be very small. The values of VLA image rms 
in a 1"7 resolution element are only 1—4 %. On the 
other hand, the absolute radio flux errors over the 
whole H II regions can be estimated from the differ- 
ences between VLA observations and single dish ob- 
servations, which are about 20 % for all sources (e.g., 
Massi et al. 1985). Finally we have added this uncer- 
tainty for that of EM for "TOTAL REGION" . The 
[Ne II] emission line photons are emitted as a result 
of collisional excitation by electrons. For all sources 
observed, the critical density {ricrit — 3.6 x 10^ cm^^ 
for the [Ne II] line) is much larger than the electron 
number density (< 10^ cm~'^) estimated from the EM 
derived from the radio flux for the angular scale of the 
MIRFI beam. Therefore we neglect the correction for 
the critical density {ne/ricrit <C 1) and the line flux 
is assumed to be proportional to EM. Then the ionic 
abundance of singly ionized neon is expressed as: 



Ne+ 



F[Ne II] (W cm-2) 



(1) 



2.3 X 10-9 rjf, Te^'"^ Q-hc/\kT^ em' 

Here F\^^^ uj is the intrinsic [Ne II] line flux corrected 
for interstellar extinction, is the solid angle of the 
observed area, and Te is the electron temperature of 
the source. 

The observed and extinction-corrected [Ne II] flux, 
and derived physical properties are listed in Table 
3. The "TOTAL REGION" means the whole region 
in which the [Ne II] emission is appreciably (3(t) de- 
tected. The "PEAK REGION" represents the 1"7 x 
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1'.'7 region (corresponding to 1 pixc;l of MIRFI) cen- 
tered at each radio continuum peak. 

Additionally we assumed electron temperature for 
optically thin H II regions, however, the extreme 
peaks of W51d and G45.12+0.13 are thought to be 
not optically thin. WC estimated peak brightness 
temperatures of 9800 K and 8100 K with correspond- 
ing optical depths of 3.71 and 1.66 for W51d and 
G45. 12+0.13, respectively (see their Table 17). They 
assiuned Tg ~ 10,000 K. Our assumptions of electron 
temperature for these two sources are 7700 K and 
8000 K. These may not be suitable for the peaks due 
to the optical depth effect, because an electron tem- 
perature distribution does make a difference from flat 
distribution around radio peak as we assumed. When 
the integrated radio fluxes before smoothing over the 
1.7 arcsec^ beam at peaks compare with that of af- 
ter smoothing, the differences are 2—12 % for each 
source. These are the additional uncertainty in their 
radio flux, succeedingly T;,, r and EM for "PEAK 
REGION". On the other hand, averages for "TO- 
TAL REGION" may not be significantly affected by 
the high optical depth at extreme radio peaks. 

We also derived the Ne+ abundance maps, us- 
ing the [Ne II] images and the radio continuum im- 
ages smoothed by the beam of MIRFI. These are 
shown in Figure 3. In order to aid the eye to see 
the trend of the Ne"*" abundance distribution, in Fig- 
ure 4 we also show the cut views through the radio 
peak. The directions of the cut are shown in Figure 3 
(A— A'). As is seen in Figure 3, the Ne"*" abundance 
has its minimum at the radio peak for all sources, 
and increases steeply towards the outer regions, espe- 
cially for W51d and G45.12-f0.13 (the ratio of the 
Nc+ abundance at the radio peak (1'.'7 x 1'.'7) to 
that of the total area ([Ne+/H+]poak/[Ne+/H+]totai) 
are 0.35 ± 0.07 for W51d and 0.28 ± 0.04 for 
G45. 12-1-0. 13), which are ionized by relatively mas- 
sive stars. On the other hand, the Ne^ distribu- 
tions are relatively uniform for G35.20— 1.74 and Mon 
R2 ([Ne+/H+]pcak/[Ne+/H+]totai are 0.87 ± 0.17 for 
G35.20-1.74 and 0.72 ± 0.12 for Mon R2), which are 
ionized by relatively low-mass stars. These rcsiilts in- 
dicate that the Ne+ distribution depends strongly on 
the spectral type of the ionizing star. 

3.3. Results for Individual Sources 

a) W51d 

WMTM reported 2.8 x 10-^* W cm-^ for the [Ne 



II] line flux observed in one raster position centered 
on the radio peak. The [Ne II] line flux (2.4 x 10"^* 
W cm~^ at the peak position) of the present work 
is in good agreement with that of WMTM within the 
flux calibration error. In this work this source was ob- 
served at many raster positions as mentioned in §2, 
and hence we can derive the detailed spatial structure 
of the [Ne II] emission and that of the Ne"*" abundance. 
Figure 2a shows that the extent of the [Ne II] emit- 
ting region is almost the same as that of the radio 
continuum emission. The observed [Ne II] line flux in- 
tegrated over the whole area is 2.8 x 10"^'' W cm~^, 
which is the highest among the sources observed. The 
cometary structure clearly visible on the radio image 
is not seen in the [Ne II] image, which may be due to 
the inferior spatial resolution of MIRFI. As shown by 
Figure 3a and Figure 4a, the Ne+ abundance at the 
radio peak is the lowest among the all sources. 
h) G45. 12+0.13 

The [Ne II] flux integrated over the whole observed 
area is 1.8 x 10~^^ W cm~^. This value is lower 
than that obtained by the previous observations with 
relatively large beams; 3.1 x 10""'^^ W cm~^ in a 7'.'5 
beam by Herter et al. (1981), and 2.6 x lO-^"^ W 
cm~^ in a 10" beam by Lester (1979). The reason for 
this discrepancy may be partly due to an insufficient 
coverage in the present work, since the [Ne II] emission 
probably extends to the outer region, especially to the 
south-west as shown in Figure 2h. So the total [Ne II] 
flux derived in this work is possibly a lower limit. 

Figure 26 shows that the general outline of the [Ne 
II] emitting region is consistent with that of the ra- 
dio continuum image, especially for the relatively low 
brightness part. Both the radio and the [Ne II] inten- 
sities generally decrease towards the north-east from 
the peaks in the south-west. Figure 3& clearly shows 
that the Ne+ abundance decreases steeply toward the 
radio peak, the same trend as that seen in W51d. 

c) GS5.20-1.74 

There have been many studies of this source as 
a part of W48 (e.g., Persi et al. 1997), however, the 
[Ne II] observation has not yet been reported. The ex- 
tended emission in the [Ne II] image is in good agree- 
ment with the 15 GHz image by WC. Both images 
show an extended emission toward the south-west di- 
rection. This suggests the ionizing sources of both the 
H II region and the Ne"*" region are the same. The to- 
tal observed [Ne II] hne flux is 9.5 x 10^^* W cm~^ 
(20.9 X 10^^^ W cm^^, after correction for interstel- 
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lar extinction). The Ne"*" abundance distribution is 
uniform over the whole area. Furthermore the Ne"*" 
abundance at the peak is larger than that of W51d 
and G45.12+0.13. 
d) Mon R2 

WC showed that the 5 GHz radio emission extends 
to more than ^ 25". The observed region by the [Ne 
II] line is a brightest part of the source (the south- 
west part; see Figure 2d). As shown in Figure 2d, the 
radio continuum and the [Ne II] emission distributions 
are in good agreement with each other. The [Ne II] 
flux integrated over the whole area is 2.0 x 10~^^ 
W cm^^ consistent with the flux of 2.3 x 10^^^ W 
cm~^ in the 7" beam obtained by Herter et al. (1982). 
The Ne"*" abundance distribution is similar to that of 
G35.20-1.74. 

4. DISCUSSION 

4.1. Effective Temperature of the Ionizing 
Star 

The Ne"*" abundance distributions differ greatly be- 
tween the sources as described in §3.2. Here we es- 
timate the effective temperature of the ionizing stars 
by comparing the observed Ne"*" abundance over the 
whole observed region with that given by H II re- 
gion models based on the latest stellar atmosphere 
models. As treated by WMTM, we use the calcula- 
tions of the line flux ratio by Stasihska and Schaerer 
(1997), which are based on CoStar stellar atmosphere 
models (Schaerer et al. 1996a, 1996b, Schaerer & 
de Koter 1997). In Figure 5 the calculated Ne"*" 
abundance is shown as a function of the stellar ef- 
fective temperature taken from grid calculations with 
Ue — 10^ cm~'^ and the LIA for elemental abun- 
dance of neon.^ From the observed Ne+ abundance 
for "TOTAL REGION" (Table 3), we estimate the 
stellar effective temperatures of the sources to be as 
follows: 39,100 tgoo (corresponding to the spectral 
type of 07.5V-08V ZAMS star) for W51d, 37,200 
tlZ K (08V-08.5V) for G45. 12-^0.13, 37,600 tf^ 
K (08V-09V) for G35.20-1.74. As for Mon R2, an 
upper limit of 33,000 K (< BOV) is derived, be- 

cause the Ne"*" abundance is as large as the LIA value 

^ The grid calculation adopted the solar abundance ([Ne/H] = 
1.7 X 10~^) as the total neon abundance (= [Ne*°'^']). We as- 
sume the calculated [Ne+/Ne*°*'''] is valid even for the case of 
the local interstellar abundance. The electron temperature cal- 
culated in this model is very close to that observed : 7500—8000 
K. 



of neon while the CoStar model grids for Toff < 33,000 
K are not available. Taking all errors into account, 
the upper limit becomes to be 34,000 K. The errors in 
the effective temperature are mostly due to the uncer- 
tainty in the absolute [Ne II] and radio fluxes. Even 
considering these errors, we can conclude that the ef- 
fective temperature of the ionizing stars in W51d is 
the highest among the sources, and that of Mon R2 
is the lowest. 

Although we estimated the effective temperature 
assuming that total abundance of Ne is LIA, the exis- 
tence of the Galactic gradient for abundances of var- 
ious ions is reported (Simpson & Rubin 1990, Henry 
& Worthey 1999). The Galactic gradient in 0/H, 
and almost certainly Nc/H as well, is —0.06 dcx/kpc. 
The range in Rq of the observed H II regions is 3.1 
kpc, which implies a difference in Nc/H of 0.18 dex 
or a factor of '^1.5. Taking these into consideration, 
total Ne abundances for each source are [Ne/H] = 
11.7 X 10"^ for W51d, 11.1 x 10"^ for G45. 12-1-0.13, 
12.5 x 10^^ for G35.20-1.74, and 8.2 x 10"^ for 
Mon R2. Using these values instead of LIA, the esti- 
mated stellar effective temperatures of the sources are 
as follows: 39,500 K (corresponding to the spectral 
type of 07.5V-08V ZAMS star) for W51d, 37,500 
K (08.5V) for G45.12+0.13, 38,200 K (< 08V) for 
G35.20-1.74, and < 33,000 K (< BOV) for Mon R2. 

The radial distribution of the Ne"*" abundance 
shown in Figures 3 and 4 can be qualitatively under- 
stood by a simple model of the ionization structure 
controlled solely by the stellar effective temperature. 
The model is described in Figure 6. In this simple 
model, the H II region is assumed to be spherically 
symmetric and the distribution of the electron num- 
ber density is assumed to be uniform over the whole 
H II region. Since the ionization potential of Ne"*" 
(41.02 eV) is higher than that of Ne° (21.56 eV), the 
ratio of amount of Ne"*" ionizing photons to that of 
Ne'^ ionizing photons depends strongly on the stellar 
effective temperature. Thus the volume ratio of the 
Ne"*""*" region to the Ne"*" region decreases as the effec- 
tive temperature decreases. The volume ratio of the 
Ne+"*" sphere to the Ne"*" shell is determined so that 
the average [Ne+"'"/Ne"'"] over the whole H II region 
is equal to that expected from the following effective 
temperature based on the model calculation shown 
in Figure 5: 39,000 K for Figure 6a, 37,000 K for 66, 
35,000 K for 6c, and 33,000 K for 6d, respectively. The 
average Ne"^ abundance distribution over the lines of 
sight crossing the equator of the sphere is also given in 
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bottom of each panel of Figure 6. The observed Ne"*" 
distributions as cut- views shown in Figure 4 are well 
explained by the model; Figure 6a corresponds well 
with the cut view of W51d, in the same way, 66 to 
G45. 12+0.13, 6c to G35.20-1.74, and 6d to Mon R2. 
However, clearly an actual H II region is not spher- 
ically symmetric, and thus the two-dimensional Ne"*" 
distribution shown in Figure 3 cannot quantitatively 
be modeled since the detailed structure of the actual 
H II region is unknown. We would emphasize, how- 
ever, the general trend seen in the variety of the Ne"*" 
abundance distributions is consistent with this simple 
model. 

An alternative explanation for the cause of the rel- 
atively low Ne+ abundance obtained for W51d and 
G45. 12-1-0. 13 may be that the [Ne II] emission is heav- 
ily extinguished by dust even though the line is in 
MIR wavelength. Especially there may be the local- 
ized variation in the extinction by dust. However, if 
the intrinsic [Ne II] fluxes from these sources are an 
order of magnitude higher that those presented in Ta- 
ble 3, then the extinction should be enormous: Ay ~ 
65 mag. Furthermore, clear near-infrared images of 
Br7 line emission which are similar to VLA images 
have been reported for both sources (W51d: Goldader 
& Wynn-Williams 1994; G45. 12+0.13: Lumsden & 
Puxley 1996). Hence it is not likely that such large 
amount of dust exists as a foreground screen. 

In Table 4, we compare the estimated spectral 
types with those derived from the FIR and the ra- 
dio luminosity. Interestingly, for the luminous sources 
W51d and G45.12+0.13, the spectral types derived 
from the FIR and the radio luminosity are much ear- 
lier than that estimated by the present work, while 
for the less luminous sources G35.20— 1.74 and Mon 
R2, both estimates are in good agreement. One pos- 
sible explanation for this trend is that the ionizing 
source is not a single star but is a star cluster, es- 
pecially for the luminous sources. The spectral type 
derived from the FIR luminosity may be earlier than 
the actual spectral type because in the relatively large 
FIR (IRAS) beam, numerous low-mass stars can con- 
tribute to the heating of the dust and hence increase 
the FIR luminosity, though they do not contribute to 
either the hydrogen or neon ionizing photons. The 
single-star spectral type derived from the radio lu- 
minosity (WC) may also be earlier than that of the 
actual one if the H II region is ionized by early type 
stars which do not emit neon ionizing photons. 

Moreover, the Ne+ abundance would not be a good 



tracer of the spectral type of the principal ionizing 
stars of the H II regions formed by very luminous 
stars like W51d and G45. 12+0.13, since Ne+ is not 
the dominant species of neon for such sources. The 
observed [Ne II] image is likely to be contaminated 
by the contribution from less luminous ionizing stars, 
leading us to imderestimate the spectral type of the 
principal ionizing star. We will further discuss this 
matter in §4.2.2, especially for G45.12+0.13. In case 
of less luminous sources, the H II regions are likely to 
be ionized by a single star, and Ne"*" is the dominant 
species of neon. Therefore, the radio-derived spectral 
type agrees well with that derived in the present work. 

The [Ne+^/Ne"*"] abundance ratio would be a bet- 
ter tracer of the spectral type of the ionizing stars, 
since this does not depend on the total elemental 
abundance of neon. However, as shown in Table 4, 
the spectral type derived by the MIR line ratio ([Ne 
II] 12.8 and [Ne III] 15.6 /xm; Simpson & Rubin 
1990) is much later than that derived by the radio lu- 
minosity. Since these MIR lines are observed with the 
large beam of IRAS LRS (25-30"), there might exist 
some contamination by less luminous ionizing stars. 
Observations of the MIR [Ne III] lines (15.6 /xm and 
36.0 iim) at higher angular resolution can provide an 
answer to this problem. 

4.2. Remarks for Individual Sources 

a) W51d 

The derived effective temperature (^ 38,500 K) is 
similar to that concluded by WMTM. Figure 3a and 
Figure 4a indicate that the Ne"*""*" is the dominant 
species of neon in W51d. 

b) G45.12+0.13 

The [Ne II] peak is not associated with the main 
radio peak (position A), but is associated with the 
second radio peak (position B; see Figure 2). Since 
the positional uncertainty of our observation is < 1", 
a spatial separation between two peaks of ^ 2" is sig- 
nificant. The Ne+ abundance derived at the [Ne II] 
peak (position B) is twice as large as that of the radio 
peak (position A). One possible explanation for this 
excess of the Ne"*" abundance at position B is that 
there is a small H II region ionized by a relatively late 
type star at position B superposed on the major H II 
region ionized by a massive star located at position A. 
Crude estimates of the [Ne II] and radio fliuces origi- 
nating from this small H II region are 2.1 x 10^^^ W 
cm~^ and 0.87 Jy, respectively, giving the Ne+ abun- 
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dance in the small H II region as 8.6 x 10^"'', which 
indicates the existence of a BO star (Teg ~ 33,000 K 
: Figure 5). However there is no firm evidence for the 
existence of the small H II region, and it is noteworthy 
that the center velocities of the [Ne II] lines at both 
positions are the same within the error of the wave- 
length calibration (< 20 km s~^). Hence the small 
H II region, if it exists, should have almost the same 
proper motion as that of the major H II region. Al- 
ternatively, the excess of the Ne+ abundance could be 
simply due to a geometrical effect. Since the effective 
temperature of the major ionizing star is relatively 
high, most of Ne"*" ions are ionized to a higher ion- 
ization state (Ne"*""*"). Then we may expect that the 
Ne+ region is distributed around the Ne++ region like 
a thin shell, and position B may be simply the tangen- 
tial component of the Ne"*" shell, enhancing the Ne"*" 
abundance. 

There are many measurements of Ne and other 
lines that can be used to model the H II regions and 
thereby determined the stellar types. Colgan et al. 
(1991) reported the results of observations of FIR 
lines. The derived effective temperature (~ 37,000 
K) is nearly consistent with our result. However 
ratio of Ne+/Ne++ is inconsistent with our results. 
The difference is caused by as follows: Colgan et al. 
used LTE models (Kurucz 1979) to estimate the ef- 
fective temperature (see Watarai et al. 1998). We 
re-estimated its using non-LTE [CoStar) models and 
the Ne+/Ne++ ratio they used, which became to be 
^ 34,000 K. Furthermore the beam size of their ob- 
servations are larger than the size of sources, hence 
there may be some possibility of contaminations. 

In order to reveal the true picture, obser- 
vations of both the [Ne II] line and [Ne III] 
line images at higher spatial resolution are essen- 
tially needed. The results will soon confirm by 
large aperture telescope observations like SOFIA 
{StratosphericObservatoryForlnfraredAstronomy) 
and SUBARU. 

c) G35.20-1.74 

Prom Figure 3c and Figure 4c, the Ne+ abundance 
is distributed uniformly like Mon R2 but is lower than 
the LIA of neon. One possible interpretation is that 
because the observed region does not cover the whole 
H II region, we may have selectively observed the ion- 
ized region in which Ne"'""'" is the dominant species 
of neon. As shown in Figure 4c, the Ne"'" abun- 
dance tends to increase gradually at the outer part. 
Furthermore, the Ne"'" abundance at the PEAK RE- 



GION is larger than those of W51d or G45.12+0.13 
(Table 3), suggesting that the thickness of Ne"'" re- 
gion toward the peak is relatively large, as shown in 
Figure 6c. Hence the Ne"'" abundance presented for 
the TOTAL REGION would be a lower limit of the 
Ne"'" abundance over the whole source and the cor- 
responding stellar effective temperature (37,000 K) 
would be an upper limit. The cut view shown in Fig- 
ure 4c corresponds well with the average Ne"'" distri- 
bution presented in Figure 6c. These results indicate 
that there exists a slightly lower temperature star (~ 
35,000 K) in G35.20-1.74, than those of W51d and 
G45.12+0.13, but higher than the one in Mon R2. 
Persi et al. (1997) proposed that these are BO— B3 
stars. WC also reported that the spectral type of the 
ionizing star is 09—09.5 based on the FIR flux and 
radio fiuxes. Further observations with much wider 
spatial coverage than the present work are necessary 
to confirm this. 
d) Mon R2 

The Ne"'" abundance map represents uniform dis- 
tribution. The value of the Ne"'" abundance is larger 
than those of the other sources and is as large as the 
LIA of neon. This implies that the spectral type of 
ionizing star is not earlier than 09. Howard, Pipher, 
& Forrest (1994) showed that IRS Isw is a B0-B3 
star, consistent with our results. Furthermore, al- 
though with a larger beam (7—30"), Herter et al. 
(1982) reported that the low excitation lines ([Ar II] 
6.99 /xm, [S III] 18.7 jiui and [Ne II] 12.8 /xm) were 
detected significantly, while high excitation one ([Ar 
III] 8.99 /im and [S IV] 10.5 /im) were not clearly de- 
tected, supporting that the spectral type is early B. 

We use Av = 23 mag as the correction for interstel- 
lar extinction (Natta et al. 1986). Natta et al. (1986) 
and Massi, Felli, & Simon (1985) pointed out that dif- 
ferent values of the extinction were derived from mea- 
surements with different beam sizes, implying that 
the extinction is not uniform. Howard, Pipher, & 
Forrest (1994) presented the extinction map from Bra 
and Br7 images, reporting an extinction ^ 27.3 mag 
at the IRS Isw- The Ne+ abundance map presented 
here shows a slight increase toward the outside from 
the peak. Hence if the value ^ 30 mag is used for 
the correction due to extinction at the peak, the Ne"'" 
abundance map would become more uniform. 
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5. SUMMARY 

Mid-Infrared [Nc II] line images of four Galactic 
UC H II regions with ^ 2" spatial resolution are pre- 
sented. The main results and conclusions are as fol- 
lows. 

1. The [Ne II] images are spatially resolved for all 
sources. In comparison with the radio VLA im- 
ages, thc^ [Sv II] linc^ and radio continuum im- 
ages are found to be generally in good agree- 
ment. We also derived the Ne+ abundance 
maps. For W51d and G45.12+0.13, the Ne+ 
abundance steeply decreases toward the radio 
peak and is lowest at the radio peak. This sug- 
gests the Ne region dominates near the cen- 
ters of these H II regions. On the other hand, 
the Ne"*" abundance is distributed rather uni- 
formly for G35.20-1.74 and Mon R2, and the 
Ne"*" abundance is almost equal to the LIA of 
neon for Mon R2. These trends indicate that 
the spectral type of the ionizing star plays a 
major role in the Ne"*" abundance distribution. 

2. We estimated the effective temperature of the 
ionizing stars from the Ne"*" abundance by us- 
ing H II region models based on the latest stellar 
atmospheric models: 39,100 tHo^ K for W51d, 
37,200 trm° ^ for G45. 12-1-0.13, 37,600 tllo° 
K for G35.20-1.74, and < 34,000 K for Mon 
R2. The range of the effective temperatures in- 
clude the effect of Galactic gradient in Ne/H. 
We further proposed a simple model of the ion- 
ization structure of neon controlled solely by 
the stellar effective temperature. The averaged 
Ne"*" distribution derived from the present ob- 
servations (Figure 4) and that from the models 
(lower graphs of each model in Figure 6) gener- 
ally agree well. As a result, the effective tem- 
perature of the ionizing star for G35.20— 1.74 is 
estimated to be ~ 35,000 K. 



can be revealed by the spatially resolved MIR line 
images. 
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3. We also presented several results and discus- 
sions for individual sources as follows; possible 
interpretations of the difference between [Ne II] 
and radio continuum peaks in G45. 12-1-0. 13, and 
the non-uniform effect of interstellar extinction 
for Mon R2. 

Finally, we have demonstrated in this paper that 
the character of the ionizing star in the UC H II region 
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Fig. 1. — Observed [Ne II] 12.814 /j,m emission spectra 
at the [Ne II] peaks of four UC H II regions, (a) W51d, 
(6) G45. 12+0.13, (c) G35.20-1.74, and (d) Mon R2. 
These spectra have not yet been corrected for inter- 
stellar extinction. It should be noted that the intrin- 
sic line fluxes are about twice as large, since Ay is ~ 
20 niag for all sources (sec text). These spectra are 
fitted with Lorentzian function (Fabry-Perot trans- 
mission profile) shown by solid line. Any continuum 
emission is not significantly observed. The error bars 
represent the la statistical error derived from the rms 
fluctuations over the observation of 40 LINE SCANs 
(512 sec in total). 

Fig. 2. — False color images of the [Ne II] emission 
intensity corrected for the interstellar extinction: (a) 
W51d, (6) G45.12-h0.13, (c) G35.20-1.74, and {d) 
Mon R2, superposed on the radio continuum images 
(contours). Radio data are taken from the 15 GHz 
observations of WC for all sources except for Mon 
R2, for which the 5 GHz observation of WC is taken. 
The peak flux density is 3.0033x10"^ Jy/beam and 
contoures are drawn at the levels —5.0, 5.0, 10, 20, 
30, 40, 50, 60, 70, 80, 90, and 95% of the peak 
for W51d; 2.4905 xlO"! Jy/beam and 1.0, 3.0, 5.0, 
10, 15, 20, 30, 40, 50, 60, 70, 80, 90, and 95% for 
G45. 12-^0.13; 8.3436x10-^ Jy/beam and 5.0, 10, 20, 
30, 40, 50, 60, 70, 80, 90, and 95% for G35.20-1.74; 
and 8.6147x10-2 Jy/beam and 2.0, 9.0, 16, 23, 30, 
37, 44, 51, 57, 64, 71, 78, 85, 92 and 99% for Mon R2. 
The observed regions in the [Ne II] line are shown as 
the colored area. The spatial sampling intervals of the 
raster mappings are typically 0'.'3 for W51d, 0'.'7 for 
G45.12-F0.13, 0'.'8 for G35.20-1.74, and 0'.'6 for Mon 
R2, respectively. 

Fig. 3. — Ne"*" abundance ([Ne^/H^]) maps of four 
UC H II regions (a) W51d, (b) G45.12-F0.13, (c) 
G35.20— 1.74, and (d) Mon R2, superposed on the ra- 
dio contours. The contour intervals are same as those 
of Figure 2. For Mon R2, the blank spot at west side 
of the radio peak is caused by the uncertainty of the 
radio contours. Thick solid lines labeled by A— A' rep- 
resent the directions of the cut views shown in Figure 
4. 

Fig. 4. Cut views of distributions of the [Ne II] and 
radio intensities as well as the Ne+ abundance. The 
directions of cuts are shown in Figure 3 (A— A'). The 
[Ne II] and radio continuum fluxes are normalized by 
the values at the peak position. Radio profiles are 



smoothed by - 2" beam (- 1 pixel of MIRFI). The 
error bars represent 1 a statistical errors of both the 
[Ne II] and the radio images. 

Fig. 5. The Nc+ abundance as a function of the 
effective temperature of the ionizing star calculated 
from the H II region model using the CoStar stellar 
atmosphere model (Stasiiiska & Schaerer 1997). To- 
tal abundance of neon (= unity at the right vertical 
axis) is assumed to be the LIA of neon (= 9.9 x 10^^ 
at the left vertical axis). The ZAMS spectral types 
corresponding to the effective temperatures are shown 
by the top horizontal axis. 

Fig. 6. The schematic model of ionization structure 
of neon using a simple and spherical H II region. The 
effective temperatures of ionizing stars correspond to 
39,000 K for (a), 37,000 K for (&), 35,000 K for (c), 
and 33,000 K for (d), respectively. The volume ra- 
tio of the Ne++ sphere to the Ne+ shell is deter- 
mined from the assiimcd effective temperature using 
the model calculation shown in Figure 5. For each 
panel, the averaged Ne+ abundance along the line of 
sight is also shown. Rc is the distance from the center 
crossing the equator projected on the sky (normalized 
by the outer radius). We assume the total abundance 
of neon is the LIA value. 
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Table 1 
TARGET LIST 



Source 


IRAS 
Designation 


a 

(1950.0) 
( h m s) 


<5 

(1950.0) 
( ° ' ") 


Distance 
from the Sun 
(kpc) 


Distance 
to the G.C. 
(kpc) 


W51d 


19213+1424 


19 21 22.29 


14 25 15.06 


7.0 


6.7 


G45.12+0.13 ... 


19111+1048 


19 11 06.23 


10 48 25.80 


9.7 


7.1 


G35.20-1.74 . . . 


18592+0108 


18 59 14.06 


01 09 03.19 


3.2 


6.2 


Mon R2 


06053-0622 


06 05 19.93 


-06 22 39.53 


0.95 


9.3 



Table 2 
OBSERVATION LOG 



Source 




Date 




Frames* 


Airmass 


Guide Star 


W51d 


Oct. 


9 


1997 


6 


1.13 - 1.89 


31 Aql. 


G45.12+0.13 


Oct. 


5 


1997 


8 


1.16 - 1.20 


7 Aql. 


G35.20-1.74 


Oct. 


5 


1997 


3 


1.49 - 1.58 


tt Aql. 




Oct. 


9 


1997 


2 


1.31 - 1.37 




Mon R2 


Oct. 


4 


1997 


3 


1.48 - 1.49 


SAO 132884 




Oct. 


8 


1997 


4 


1.48 - 1.58 






Oct. 


9 


1997 


8 


1.48 - 1.78 






Oct. 


10 


1997 


6 


1.48 - 1.67 





* See text for the definition. 
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Table 3 

DERIVED PHYSICAL PROPERTIES 



F[Ne II] flux Radio parameters Abundance 



Source Integ. Area observed corrected " Radio fiux T;, t EM [Ne+/H+] 

(10"" sr) (10-^** W cm"^) (mjy) (K) (lO'' pc cm-'') (10"^) 

TOTAL HECIUN ^ ~ 

W51d 164 27.9 ± 6.6 55.3 ± 13.0 13900 ± 1400 1240 ± 120 0.18 ± 0.02 11.4 ± 1.1 1.3 ± 0.3 

G45.12+0.13 155 18.2 ± 4.2 42.7 ± 9.8 5210 ± 520 490 ± 50 0.06 ± 0.01 4.3 ± 0.4 2.8 ± 0.7 

G35.20-1.74 95 9.5 ± 2.3 20.9 ± 5.1 3140 ± 310 480 ± 50 0.07 ± 0.01 4.1 ± 0.4 2.3 ± 0.6 

Mon R2 257 20.3 ± 4.7 44.5 ± 10.3 1610 ± 160 870 ± 90 0.12 ± 0.01 0.7 ± 0.1 10.7 ± 2.7 

PEAK REGION ^ 

W51d 6.8 2.4 ± 0.6 4.7 ± 1.1 2340 ± 210 5040 ± 470 1.06 ± 0.18 66.8 ± 11.2 0.45 ± 0.13 

G45.12+0.13 A . . . 6.8 1.2 ± 0.3 2.7 ± 0.6 1050 ± 130 2250 ± 270 0.33 ± 0.05 22.5 ± 3.2 0.79 ± 0.21 

B... 6.8 1.5 ± 0.3 3.5 ± 0.8 750 ± 50 1610 ± 120 0.23 ± 0.02 15.3 ± 1.3 1.5 ±0.4 

G35.20-1.74 6.8 1.4 ± 0.3 3.1 ± 0.8 520 ± 60 1130 ± 130 0.16 ± 0.02 10.2 ± 1.3 2.0 ± 0.5 

Mon R2 6.8 1.0 ± 0.2 2.2 ± 0.5 100 ± 2 2060 ± 40 0.32 ± 0.01 1.9 ± 0.04 7.7 ± 1.8 



Note.— We use the electron temperature as follows; W51d : 7,700K (Mehringer 1994), G45. 12+0.13 : 8,000K (Lumsden & Puxley 1996), G35.20-1.74 : 
7,500 K (Wink, Wilson, & Bieging 1983), Mon R2 : 7,600 K (Dowens et al. 1975). 

^ Corrected flux for interstellar extinction. 

^ Radio continuum fluxes at 15 GHz for all sources except for Mon R2. As for Mon R2, 5 GHz flux is adopted (Wood & Churchwell 1989). 

■= TOTAL REGIONS are defined as the regions where the [Ne II] flux is appreciably detected (above 3<t). PEAK REGIONS correspond to the MIRFI 1 
pixel regions (1'.'7 X 1'.'7) at the radio peaks except for G45.12+0.13 B, which is the [Ne II] emission peak. 



Table 4 

SPECTRAL TYPE OF IONIZING STAR 



Source FIR Luminosity ^ Radio Luminosity ^ MIR Line Ratio This Work 

Single Star Cluster Single Star 

W51d 04 04.5 05.5 07.5V-08V 

G45.12+0.13. . . 04 05 06.5 09.5V 08V-08.5V 

G35.20-1.74. . . 06 09.5 09 08V-09V 

Mon R2 BO BOV <BOV " 



IRAS PSC (Wood & Churchwell 1989). 

^ Converted from A^Lyc luminosity (Wood & Churchwell 1989). 

Converted from [Nc+/Ne++] given by Simpson & Rubin 1990 using the calculations based on the CoStar 
stellar atmosphere model (Stasinska & Schaerer 1997). 

Downes et al. 1975. 
" The CoStar model grids for Tefi < 33,000K (< BO) are not available. 



13 




Wavelength 





Right Ascension (18 h 59 i^i) 



Riglit Ascension (06 h 05 m) 
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Spectral Type 
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